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Aims. The aim of this work is to investigate in a physical and quantitative way the spectral evolution of bright Neutron Star Low-Mass X-ray 
Binaries (NS LMXBs), with special regard to the transient hard X-ray tails. 

Methods. We analyzed INTEGRAL data for five sources (GX 5-1, GX 349+2, GX 13+1, GX 3+1, GX 9+1) and built broad-band X-ray 
^— «j spectra from JEM-X1 and IBIS/ISGRI data. For each source, X-ray spectra from different states were fitted with the recently proposed model 
£ — , compTB. 

£SJ ' Results. The spectra have been fit with a two-compTB model. In all cases the first compTB describes the dominant part of the spectrum 
that we interpret as thermal Comptonization of soft seed photons (<1 keV), likely from the accretion disk, by a 3-5 keV corona. In all cases, 
this component does not evolve much in terms of Comptonization efficiency, with the system converging to thermal equilibrium for increasing 
accretion rate. The second compTB varies more dramatically spanning from bulk plus thermal Comptonization of blackbody seed photons to 
, the blackbody emission alone. These seed photons (R< 12 km, kT s > 1 keV), likely from the neutron star and the innermost part of the system, 
• • ■ the Transition Layer, are Comptonized by matter in a converging flow. The presence and nature of this second compTB component (be it a 
, ^ pure blackbody or Comptonized) are related to the inner local accretion rate which can influence the transient behaviour of the hard tail: high 
K> values of accretion rates correspond to an efficient Bulk Comptonization process (bulk parameter 8^0) while even higher values of accretion 
» I ■ rates suppress the Comptonization, resulting in simple blackbody emission (S — 0). 

03 Conclusions. The spectral evolution of the sources has been successfully studied in terms of thermal and Bulk Comptonization efficiency, in 
relation to the physical conditions in the Transition Layer. 

Key words. Stars: individual: GX 5-1, GX 13+1, GX 3+1, GX 9+1, GX 349+2 - X-rays: binaries - binaries: close - stars: neutron - 
accretion 



1. Introduction 

Low-Mass X-ray Binaries (LMXBs) are systems where a com- 
pact object, either a neutron star (NS) or a black hole candi- 
date (BHC), accretes matter via Roche lobe overflow from a 
normal companion of mass M < 1M Q ; a peculiar characteris- 
tic to this type of system is the formation of an accretion disk 
in the orbital plane near the compact object. In this paper we 
studied persistently bright NS LMXBs (L x ~ 10 37 - 10 38 
erg sec -1 ) granting long-lasting observability with instruments 
operating in the soft/hard X-ray range. The spectra of these 
sources are usually described as the sum of two components 



(e.g lMitsuda et alll984IWhite et alJll988l:lBarretll200ll) : a soft 
component often likely associated with the accretion disk or NS 
and a hard component interpreted as thermal Comptonization 
of soft seed photons from the disk and/or the NS by high 
temperature plasma of electrons near the compact object (so- 
called corona). The advent of broad-band X-ray missions, such 
as BeppoSAX, RXTE, INTEGRAL, revealed the presence of 
a spectral hardening (so-called "hard tails") above ~30keV 



on top of otherw ise soft spectra (e.g lDi Salvo & Ste lla 2002; 
Paizi s et al . 2006, hereafter P06, and references therein). These 



hard tails, mostly fit with phenomenological m odels such as a 
power-law, have been detected in Z sources (e.g. lFrontera et alj 
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19981: iDi Salvo et al.ll2002l) and also in the bright Atoll source 
GX 13+1 (P06). 

To explain the origin of these hard tails, different models 
have been proposed across the years, such as direct synchrotron 
emission from a jet ( Markoff et al. 2005b. hybrid thermal/n on- 
thermal Comptonisation dCoppill 1 9991: IDi Salvo et al.ll2006l) or 
more recently, bulk motio n Comptonisation (ITitarchuk et al.l 
1 19971 iFarinelli et al.ll2008l hereafter TMK97 and F08, respec- 
tively). The adoption of different models which satisfactorily 
fit the hard X-ray component in NS LMXBs clearly shows that 
it is necessary to look for other observable quantities in these 
systems, such as their timing properties or correlation among 
the spectral parameters. 

For black-hole sources, a fundamental step forward came 
with the unambiguous discovery of the saturation of the spec- 
tral index (as a function of either the low QPO frequency or 
seed photon BB normalization) w hen sources moves from the 
low/h ard to the high/soft state fe.g. lShaposhnikov & Titarchuk 
20091 and references therein). Such a saturation can be natu- 



rally explained in terms of the presence of a converging flow 
close to the compact object. 

For NS systems, the situation is more complex because of 
the presence of a firm surface in the compact object that plays a 
fundamental role in determining the hydrodynamical and radia- 
tive system configuration with, in turn, a different spectral ap- 
pearance. However, given that a PL-like emission can be seen 
in both BH and NS systems and on the basis of the above re- 
ported results for BHs, it is a reasonable working hypothesis 
that converging flow (with unavoidable qualitative and quan- 
titative differences) is a characteristic of NS systems too, at 
least in some particular accretion rate states. Close to the NS 
surface, the radial component of the velocity of matter (Vr) 
dominates over the azimuthal one and it is straightforward to 
find that Vr ~0.56 c at Rns, factor that is actually reduced 
by the presence of pressure gradients (due to radiation, gas and 
magnetic field). Our interpretation is that multiple thermal plus 
bulk Compton scatterings, that occur in a high optical depth 
environment, are at the origin of the hard tails observed in NS 
LMXBs, extending over the thermal continuum up to at least 
lOOkeV. A self-consistent physical treatment of the innermost 
region in the case of neutron stars is presently being devel- 
oped by Titarchuk & Farinelli (in preparation, see Sect. 14. U . In 
this paper we adopt the bulk motion Comptonization scenario 
which also enabl es a homogeneous interpre t ation with previous 
works (e.g. P06. IFarinelli et al.ll2007[l2008ll2009h . 

P06 proposed for the first time a unified physical sce- 
nario to explain the spectral evolution of NS LMXBs, includ- 
ing the peculiar transient hard tail. All the observed spectral 
states could be well fit in terms of the interplay of thermal 
and bulk Com ptonization (TC and BC, resp ectively) using the 



BMC model ffitarchuketai]ll996l 1 19971 hereafter TMK96 



and TMK97) in XSPEC. Thermal Comptonization operates in 
a region where a hot plasma is present and acts on a population 
of soft seed photons, hardening the injected seed photon spec- 
trum (from the accretion disk, Fig . |T|) ; the resulting spectrum is 
extended to higher energies depending on the Comptonization 
efficiency. Bulk Comptonization operates in the inner part of 
the system between the Keplerian accretion disk and the NS 
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Fig. 1. Schematic view of the proposed scenario for thermal 
and bulk Comptonization regions in LMXBs hosting NSs. 

surface, the Transition Layer, hereafter TL (Fig. Q]). In the in- 
ner part of the TL, matter falls towards the NS surface with a 
preferred (bulk) motion, acting on the inner seed photons of the 
system; the result is a power-law-shape spectrum at high ener- 
gies due to the energy transfer from the matter kinetic energy 
to the photon field (TMK96, TMK97). P06 found evidence 
that bulk motion is always present and its strength is related 
to the accretion rate and is suppressed only in the presence of 
high local luminosit y. This scenario has been confirmed also by 
(IFarinelli et al.l 120071 hereafter F07), who re-analyzed broad- 
band BeppoSAX data, applying the BMC model to GX 17+2. 

In the quest to study in a quantitative way the evolution of 
the parameter s describing the inne rmost physical conditions of 
NS LMXBs, IFarinelli et all d2008l) (hereafter F08), developed 
a new Comptonization model, compTB^. We recall here the 
basics of this recent model, referring the reader to F08 for a 
more detailed description. 

The total emerging spectrum is given by 



F(E) 



Cn 



{BB 

(a) 



A x BB * G) 

(b) 



(1) 



where C n is the normalization constant, A is the illumination 
factor, C n BB/(A + 1) is the seed photon spectrum directly 
seen by the observer and not modified by Comptonization pro- 
cesses. C n Aj (A + 1) x BB * G is the Comptonized spectrum 
obtained by convolution of a seed photon spectrum (blackbody, 
BB) with the Green's Function^, in order to evaluate the effect 
of TC and BC on the seed photon field. 

Parameters of the compTB XSPEC model are the seed pho- 
ton temperature kT s , the electron plasma temperature kT e , the 
spectrum energy slope (i.e. overall Comptonization efficiency) 
a (photon index T = a + 1), the bulk parameter <5 that quan- 
tifies the efficiency of BC over TC and log A, which assigns a 
different weight to the two components (a) and (b). This model 
enables the co-existence of the direct seed photon component 
and its Comptonized part, all obtained in a self-consistent way. 
In Eq.[Tl we note that for log A = —8, we have only (a), i.e. 
the direct seed photon component, while for log A = 8, the di- 
rect component is no longer visible and we have only thermal 



1 http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/models/comptb.html 

2 The Green's Function, G, is the response of the system to the in- 
jection of a monochromatic line, see F08 for details. 
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Table 1. The five Bright NS LMXBs studied in t h is pa- 
per. References for di stance and N H : Grimm et all (2002 ) 
and lUeda et alj d2005l) (GX 5 -1). iGrimmetalJ d2002l) and 
Iaria etalJ d20C^h (GX 349+2), iBandvopadhvav et al] d 19991) 
and lUeda et all d2004l) (GX 13+ 1), iKuulkers & van der KlisI 
(20001) and lChenevez et al.l (Eo06) (GX 3+ll. llaria et alJd2005l 



(GX9+1). 



Source 




K 


pointings # 


Z sources 


GX5-1 


7.2 


2.8 


137 


GX 349+2 


9.2 


0.77 


51 


Atoll sources 


GX 13+1 


7 


3.2 


35 


GX3+1 


5 


1.66 


166 


GX9+1 


5 


0.8 


87 



a The distance is in kpc. 
b In units of 10 22 cm" 2 . 



2. Observations and data analysis 

2. 1. Analysis process 

The sample of the NS LMXB chosen is given in Tab. [T] and 
comprises five bright ("GX") bulge sources. We have ana- 
lyzed INTEGRAL pointings collected from September 2006 
to September 2007 from the A04 and A05 INTEGRAL Key 
Programmes. The exposure time for each pointing is variable 
(At ~ 2000 sec). The selected pointings satisfy the condition 
that every source is < 3° from the pointing direct ion. Analysis 



was performed on JEM-X1 and IBIS/ISGRI dLebrun et al 



2003I) data with the software package OS A, version 7.0. 

The energy bands selected for JEM-X1 imaging analysis 
are: 3-10keV and 10-20keV while the energy band selected 
for IBIS/ISGRI imaging analysis is: 20^40 keV. The mosaic of 
the Galactic Centre (Fig. [2j was constructed using all available 
pointings in the IBIS/ISGRI energy band: as expected, all the 
sources studied were active and detectable. We used the official 
INTEGRAL catalog to extract spectra for the single sources and 
the JEM-X default response matrix (256 bins) and a new re- 
binned (57-bins) IBIS/ISGRI response matrix. 



2.2. Spectral classification with HR 



plus bulk Comptonization. In the case of 5 = 0, the BC contri- 
bution is neglected, retai ning only the TC effects (equivalent to 
compTT, iTitarchukll 1 994 . 

The first application of the compTB model and of the pro- 
posed physical scenario to the spectral evolution of a single 
NS LMXB, to trace in a quantitative w ay the evolution of the 



physical parameters, was presented in iFarinelli et afl d2009h 



(hereafter F09) on Cyg X-2, using BeppoSAX data. F09 fit- 
ted the spectra of the source with two compTB models. The 
first one, with bulk parameter (5=0, was the dominant com- 
ponent and described a pure thermal Comptonization process. 
This component, due to the Comptonization of cold disk pho- 
tons by the outer TL, is roughly constant with time and the 
accretion rate. The second compTB model describes the over- 
all Comptonization (thermal plus bulk, variable 8 parameter) of 
hotter seed photons close to the neutron star surface by the in- 
ner region of the TL (Fig.Q]). This component is highly variable 
with time (and local accretion rate) ranging from significant 
Comptonization of the hot seed photons (thermal and bulk) to 
a simple blackbody-like spectrum. 

In this work, we extend the above study to five persistently 
bright NS LMXBs (GX 5-1, GX 349+2, GX 13+1, GX 3+1 
and GX 9+1). The choice of sources was driven by the valid- 
ity range of the compTB model (diffusion regime, see F08), 
associated to a high accretion rate where the transient hard 
tail can be appreciated. The selected sources are located in the 
Galactic Centre where, because of the presence of a large num- 
ber of X-nr^sources, ^ngular resolution is an issue. The imager 
IBIS dUbertiniet alj|2O03l) and the X-ra y instrument JEM-X 
dLundet alj|2003h on-board INTEGRAL dWinkleret alj|2003l) 
have the necessary angular resolution to disentangle the differ- 
ent contributions, allowing to build good quality broad-band 
X-ray spectra (3-200 keV). 



In order to study the spectral evolution of the sources, we as- 
sociated a numeric parameter to the spectrum of each pointing, 
the Hardness Ratio (hereafter HR). The HR is computed from 
the counts per second extracted from the two energy bands of 
JEM-X analysis as follows: 



(2) 



where S (soft) is the count rate for the 3-10keV energy band 
while H (hard) is the count rate for the 10-20keV energy band. 
HR is an indicator of the source spectral state: the lower the HR 
the harder the spectrum. Despite its unusual/counterintuitive 
form, we adopted the definition (O because the good statis- 
tics showed by JEM-X granted us integer values for a more 
hand-light management of HR boundaries and, as explained in 
section |4~T1 we wanted to connect the increase of HR with the 
increase of the accretion rate. 

Considering the overall behaviour, each source presents a 
characteristic HR variability spanning wide ranges (for exam- 
ple see Fig. [3] for GX 5-1). From a morphological evaluation 
of the spectral characteristics together with HR values, we clas- 
sified every spectrum and averaged similar spectra to improve 
the overall statistics. 

2.3. The spectral fitting 

Before spectral fitting, all the spectra were re-binned so as to 
contain at least 25 counts in every channel; we considered only 
the counts within the interval 5-18 keV for JEM-X1 spectra 
and the counts starting at 22 keV for IBIS/ISGRI. Systematic 
errors of 2% for IBIS/ISGRI data and 3% for JEM-X1 data 
were added. Instrumental cross-calibration factors were al- 
lowed to vary between 0.8 and 1.2 (IBIS/ISGRI) with JEM-X1 
frozen at 1 . In all cases the value of the constant stabilized itself 
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Fig. 2. IBIS/ISGRI 20-40 keV mosaic in Galactic coordinates of the Galactic Bulge (about 3.3 Msec). Only the sources studied 
in this paper are labeled. 
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Fig. 3. Example of the HR variation for the source GX 5-1 con- 
sidering all available pointings. 

at the maximum value, 1 .2. For the spectral analysis, we used 
XSPEC version 12. 

In all cases we started the fitting process with the simplest 
general model considering only one compTB and then in more 
complex cases, we added a second compTB, as in F08. Hence 
(excluding the cross-calibration constant), in the most complex 
cases, we used: 

F{E) = wabs * (compTB[l] + compTB[2]). (3) 

where wabs dMorrison & McCammon|[l983 ) is the multiplica- 
tive model that takes into account the interstellar absorption at 



low X-ray energies; the column-density, Nh , was fixed during 
the fitting to the values shown in TableQ] as the absence of data 
below 5 keV does not allow the correct determination of this 
parameter through the fitting process. 

We note that compTB[l], with parameters set to logA = 
8 and S = 0, is a thermal compTB accounting only for TC 
processes (equivalent to compTT). Instead, compTB[2], when 
present, is either a mixed compTB (shape-wise equivalent to « 
bbody + PL, where we consider both TC and BC processes), or 
a bbody (neither TC nor BC setting log A = —8). For the cases 
in which the fit required two compTB components, a single 
temperature of the plasma kT e was considered (see section [3] 
for details). 

2.4. Other parameters 

Some other quantities can be obtained from the fit processes 
in order to have more information about the system and the 
environment around the NS: 

- Ftot: the unabsorbed total flux (erg sec -1 cm~ 2 ) in 
the interval 0. 1-200 keV and relative luminosity Ltot — 
AirD 2 Ftot (using distances in Table [T); 

- CAF: the Compton Amplification Factor that allows us to 
evaluate the total energy gain of the seed photons through 
the TC process. It is defined as CAF = fluxl/flux2, 
where fluxl (erg sec -1 cm~ 2 ) is the unabsorbed energy flux 
associated to the thermally Comptonized photons (calcu- 
lated on compTB[1] with Nh = 0), and flux2 (erg sec -1 
cm -2 ) is the unabsorbed energy flux associated to the seed 
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photons (calculated on C0impTB[1] with log A = —8 and 
N H = 0) in the energy interval 0. 1-200 keV; 
The optical depth r: the Comptonization parameter a from 
the fit is linked to kT e and r as: 



Table 2. The number of pointings and exposition time for each 
HR interval for the NS LMXBs studied in this paper. 




C T fcT e (r+|) 2 



(4) 



where C T — 3 for a spherical geometry with opti- 
cal radius r, while C T = 12 for a "slab" geometry 
(Titarchuk & Lvubarskiil[l995l) . n this paper we presumed 
a slab geometry corresponding to a not very extended 
corona (see also F08); 

Rbt>: the blackbody radius of compTB[2], computed using 

L bb = AnR 2 bb aT A . 



3. Results 

The detailed results per source are given in the next subsections 
while here we note a few general trends and issues. Each source 
shows a given range of variability for HR values which enabled 
us to identify a different number of intervals (see Table 13 and 
spectral states per source . Our sources showed all types of be- 
haviours, and hence we used a model with: two compTBs in 
which the first is thermal and the second mixed (GX 5-1 and 
GX 13+1); two compTBs in which the first is thermal and the 
second a blackbody (GX 9+1 and GX 349+2); one thermal 
compTB (GX 3+1). The absence of data below 5keV led to 
a high uncertainty in determining kT s in the compTB[l] com- 
ponent (the softest). In fact, for the majority of cases only an 
upper limit could be determined. 

As already experienced in F08 and F09, in the cases where 
we applied a mixed compTB[2], the lack of an observable cut- 
off prevented us to allow a, 8 and log A to be free during 
the fit: indeed if the number of the free parameters exceeds 
the number of the observable quantities (normalization, cut- 
off and the PL slope), a degeneration of model parameters oc- 
curs and hence more than one model solution can be obtained. 
The degeneration can be reduced: as in F08, the kT e value 
for compTB[2] was set equal to the compTB[l] one, consider- 
ing a single plasma temperature (indeed a preliminary fit pro- 
vided similar values for both kT e ). Furthermore, we fixed log A 
(GX 5-1) or a (GX 13+1) to the values corresponding to a x 2 
minimum (chosen inside an interval calculated with the XSPEC 
command steppar). We left 8 free during the fit to evaluate the 
bulk efficiency. 



3.1. GX5-1 

The variability range of HR (definition 0) spans between ~5- 
18 (as shown in Fig. [3). By classifying the spectral shapes, 
we obtained four mean spectra : the first spectrum, the hard- 
est one, corresponding to HR<7, shows, as expected, a hard 
tail between 40 and 150keV. Indeed a hard tail in the spectra 
of GX 5-1 had already been detected by previous observations 
dAsai et alJll994 IPaizis et al.ll2005i P06). The second spec- 



Interval 


pointings 


T 


(ksec) 




N 


JEM-X 


ISGRI 


GX5-1 


HR<7 


36 


102 


69 


7<HR<10 


52 


155 


106 


10<HR<12 


23 


58 


38 


HR>12 


26 


80 


53 


GX 13+1 


HR<10 


11 


29 


18 


10<HR<13 


10 


26 


17 


HR>13 


14 


35 


23 


GX9+1 


HR<7.5 


19 


50 


34 


7.5<HR<8.5 


47 


119 


79 


HR>8.5 


21 


54 


37 


GX 349+2 


HR<7 


20 


52 


35 


HR>7 


31 


82 


53 


GX3+1 


HR<8 


37 


107 


70 


8<HR<9.5 


103 


323 


210 


HR>9.5 


26 


78 


52 



trum (7<HR<10) shows a weaker hard tail between 40 and 



150 keV; the third one (10<HR<12) does not show a hard tail 
and is soft; the last s pectrum is the softest (HR>12), again with 
no hard tail (see also lKuulkers et a l. 2007; Paizis et al.l2005l on 
the variability of GX 5-1). 

As required by the data, we used the two-component model 
for all spectra in the fitting process and obtained the following 
results (see Table[3]l: 

- with the compTB[l] component, we describe the presence 
of the TC process as the dominant mechanism in all the 
four spectra. For higher values of HR, we observe that the 
CAF decreases, while a increases. This is consistent with a 
reduction in the efficiency of the TC process; 

- the seed photon temperatures kT s of the two components 
are different: in the thermal compTB[l], kT s remains al- 
ways below IkeV while in the mixed compTB[2] it has 
higher values > 1.5 keV; 

- considering the errors associated, the plasma temperature 
kT e remains within the interval 3^1 keV; 

- with compTB[2] in all spectra, we describe an additional 
component that shows a specific evolution: in the first two 
spectra (HR<7, 7<HR<10) we obtained a mixed compTB 
that allows us to evaluate the BC contribution with the pa- 
rameter 8 ^ 0. The spectra are representative of the inter- 
mediate state (according to the definition given by P06). For 
increasing values of HR, log A shows a decreasing trend 
that assigns less importance to the mixed component. In 
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Table 3. Best-fit parameters of the multi-component model 
wabs • (compTB[l] + compTB[2]) for the source GX 5—1. 
Errors are computed at 90% confidence level for a single pa- 
rameter. 



Table 4. Best-fit parameters of the multi-component model 
wabs - (compTB[l] +compTB[2]) for the source GX 13+1. 
Errors are computed at 90% confidence level for a single pa- 
rameter. 













parameter 


HR<10 


10<HR<13 


HR>13 


parameter 


HR<7 


7<HR<10 


10<HR<12 


HR>12 




compTB[l] (thermal: log A = 8, S 


= 0) 




compTB[l] (thermal: log A 


= 8, 5 = 0) 




kTs a 


0.2 (< 1.8) 


u - yz 0.08 


0.893^' 


kT s a 
kT e a 
a 


0.4 (< 0.6) 

o 4 t-+0.05 
1 C i+ ' :L 


0.6 (< 0.7) 

q 11+0.05 
°' ±± -0.13 

1 S +0 1 
J "°-0.2 


0.4 (< 0.6) 
3.7±g;g 

o -i +0.2 

2 - 1 Io.5 


* +1) 1 

0.4^i 

[4] 
2 -5-o.i 


kT e ° 
a 

d 

T 


O+0.4 
■J-0.3 

3.8 


+ 12 
•J-5-0.6 
9 4 +0.2 
Z -*-0.4 

2.4 


4-211;? 

q-t-u.o 
°-0.5 

1.7 


r d 


3.6 


3.8 


2.6 


2.1 


CAF d 


1.8 


1.2 


1.1 


CAF d 


1.6 


1.5 


1.4 


1.3 




com 


pTB[2] 








compTB[2] 






kTs a 


1.7 (> 0.7) 






kT s a 
R bb M 
kT e a 
a 


[kT e ] 

1 9+0-9 

J -*-o.a 


1 6 +ua 

D - z -0.8 

[kT e ] 

1 1+ - 4 

L - L -0.3 


1 71+0.U4 
± - ' 1 -0.04 

5.6±g;i 


1 e 7 +0.U2 
1 -°'-0.02 


P b,d 

kT e a 
a 
5 


2(<H) 
[kT e ] 
[0.78] 
103 (> 12) 






S 


36 (> 18) 


33 (> 14) 






log A 


-1.3 (> -1.9) 






log A 


[-1.82] 


[-2.19] 


[-8] 


[-8] 


T c ;d 

J^tot 


1.8 


0.6 


0.5 


r c.d 
J^tot 


2.5 


2.5 


3.5 


3.5 


T c.d 
J^bb 


0.03 






T c >d 
Lbb 


0.3 


0.3 


0.4 


0.5 




(1.7%) 








(12%) 


(12%) 


(11%) 


(14%) 




14/20 


32/19 


22/19 


X Z /dof 


133/130 


101/118 


82/114 


96/109 











a In keV. 

b In km. 

c In units of 10 38 erg sec -1 , in the 0. 1-200 keV energy range. 

d Computed as reported in section [2~4l 



a In keV. 

b In km. 

c In units of 10 38 erg sec -1 , in the 0.1-200 keV energy range. 

d Computed as reported in section |2l4l 



fact, in the third and fourth spectra (10<HR<12, HR>12), 
we found a best-fit using the compTB[2] in the blackbody 
mode (log A = —8 ); 
- with increasing values of HR, the size of the blackbody- 
emitting region, Rhb in compTB[2], increases, remaining 
compatible with the NS and/or TL dimensions. 

In Fig. |4] we show the EF(E) spectra (keVcm -2 sec -1 ) of 
GX 5-1 for two of the four selected HR intervals (correspond- 
ing to the hardest and the softest case). 

3.2. GX 13+ 1 

Examining the shape of the GX 13+1 JEM-X1 and 
IBIS/ISGRI spectra, we identified three intervals in which we 
recognized different phases of the spectral evolution and so we 
obtained three mean spectra : the first mean spectrum corre- 
sponding to HR< 10 shows a hard tail between 30 and 150keV 



cnaractenzing an intermediate state, /v nara tan in tne spec- 
tra of GX 13+1 had been discovered using INTEGRAL data 
by P06. The mean spectrum corresponding to 10<HR<13 
shows a high/soft state, deprived of the hard tail, similar to 
the mean spectrum corresponding to HR>13. Spectral analysis 
was carried out using a double-component model only for the 
HR<10 spectrum, while the other two were fit with one ther- 
mal compTB (see Table|4] Fig. [5J. Looking at the parameters 
we note: 



- the compTB[l] component, present in all spectra, describes 
the TC process. Again, for increasing values of HR we ob- 
serve a strong decrease of CAF values and increase of a, a 
signature of decreasing efficiency of the TC process; 

- the seed photon temperatures kT s in the two components 
are different in the first spectrum (in the other two spectra 
only one thermal compTB was needed); 

- the plasma temperature kT e remains within the range of 3- 
4keV; 

- we applied the compTB[2] component only to the first 
spectrum as a mixed compTB: we detected a hard tail and 
an active BC process with 5^0; the a value was fixed to 
a value that corresponds to a minimum x 2 ( see section[3]i; 

- the luminosity of the blackbody Lhh is very low with re- 
spect to the total X-ray luminosity of the source; the di- 
mension of the blackbody emission region, Rhb ~2 km, is 
compatible with part of the NS surface and/or TL dimen- 
sions. 



3.3. GX9+1 

We classified every JEM-X1 and IBIS/ISGRI spectrum of 
GX 9+1 according to their shape and identified three HR in- 
tervals (HR<7.5, 7.5<HR<8.5 and HR>8.5) obtaining three 
mean spectra. All spectra showed the high/soft state, in fact, 
the spectral analysis was carried out with a double-component 
model using a thermal compTB plus a second compTB as a 
pure blackbody (see Table|5] Fig. |6j: 
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Fig. 4. Unabsorbed EF(E) spectra of GX 5-1 for the HR intervals: HR<7 (left), and HR>12 (right). The dotted lines are the two 
separated compTB components: dash for compTB[l], dash-dot compTB[2]. 
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Fig. 5. Unabsorbed EF(E) spectra of GX 13+1 for the HR intervals: HR<10 (left), HR>13 (right). The dotted lines are the two 
separated compTB components: dash for compTB[l], dash-dot compTB[2]. 



- the TC process was described with compTB[l] in all three 
spectra; 

- similarly to the previous cases, the seed photon tempera- 
tures kT s differ in the two components: about 0.6 keV for 
the first one and near 2 keV for the second one; 

- in this case, the plasma temperature kT e is not well con- 
strained by the data; 

- the presence of the second component compTB [2] re- 
veals a constant blackbody contribution. Below 20 ke V the 
blackbody emission is more important than the TC one. 
The size of the blackbody-emitting region is nearly con- 
stant (~3.4 km), as is the corresponding X-ray luminosity 
(0.2 x 10 38 erg sec -1 ). Despite the different HR selection, 
these spectra do not show a dramatic evolution, as shown 
in the previous cases. 



3.4. GX 349+2 

For GX 349+2, we identified two intervals (HR<7 and HR>7) 
showing similar spectra from which we obtained two mean 
spectra: both are of the high/soft state well described with a 
double-component model consisting of a thermal compTB and 
a second compTB as blackbody. We obtained the following re- 
sults (see Table [6] Fig.|7]i: 

- with compTB[l], we describe the TC process; for higher 
values of HR we do not observe the typical a growth but 



Table 5. Best-fit parameters of the multi-component model 
wabs • (compTB[l] + compTB[2]) for the source GX 9+1. 
Errors are computed at 90% confidence level for a single pa- 
rameter. 



parameter 


HR<7.5 


7.5<HR<8.5 


HR>8.5 




compTB[l] (thermal: log A = 8,5 


= 0) 


kT s a 
kT e a 
a 


0.3 (< 0.6) 

4 +3 

9 1+0.6 
^•'--0.6 


0.3 (< 0.6) 

5+i 2 

9 4+0-7 


0.3 (< 0.6) 

5(>3) 
2 4+ - 9 




2.5 


2.0 


1.9 


CAF d 


1.5 


1.4 


1.4 




compTB[2] (blackbody: log A = 


-8) 


kT 3 a 
R bb M 


1 - ay: -0.04 

3.5l° :l 


1 SQ+U.03 

O q + 0.1 

°-°-0.1 


J-^-o.os 
3.4i° ;* 


kT e a 








a 








8 








log A 


[-8] 


[-8] 


[-8] 


T c ;d 

J^tot 


1.1 


1.3 


1.1 


T c ,d 
J^bb 


0.2 


0.2 


0.2 




(18%) 


(15%) 


(18%) 


X'fdof 


67/106 


83/103 


79/105 



a In keV. 

b In km. 

c In units of 10 38 erg sec -1 , in the 0.1-200 keV energy range. 

d Computed as reported in section |2l4l 
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Fig. 6. Unabsorbed EF(E) spectrum for the HR interval 
7.5<HR<8.5 for GX 9+1. The dotted lines are the two sep- 
arated compTB components: dash for compTB[l], dash-dot 
compTB[2]. 

Table 6. Best-fit parameters of the multi-component model 
wabs • (compTB[l] + compTB[2]) for the source GX 349+2. 
Errors are computed at 90% confidence level for a single pa- 
rameter. 



parameter 


HR<7 


HR>7 




compTB[l] (thermahlog 


A = 8, 8 = 0) 


kTs 
kT e 
a 


0.3 (< 0.9) 

q+1.8 
°-0.1 
1 I" 1 " 02 

5.0 


0.4(< 2.2) 

9 S+ 4 ' 7 
^•°-0.08 
-, -,+0.5 

5.2 


CAF d 


2.1 


2.0 


compTB[2] (blackbody: log ,4 = -8) 


kTs 
R bb b - d 


1 y+0.3 
1 - > -0.3 

5.2±?:« 


1 2 +u - Y 


kT e 






a 






S 






log A 


[-8] 


[-8] 


T c,d 
■utot 


3.4 


2.2 


J^bb 


0.3 


0.3 




(9%) 


(14%) 


x'/dof 


104/115 


76/114 



° In keV. 
6 In km. 



In units of 10 38 erg sec 1 , in the 0. 1-200 keV energy range. 
Computed as reported in section l2"l4l 



instead it is constant, which brings to a constant CAF value 
and invariant TC efficiency; 

similarly to previous sources, the seed photon temperatures 
are different when two compTB components are taken into 
consideration: kT s < IkeV for the thermal compTB, 
higher temperature for the blackbody photons; 
the plasma temperature stays always around 3 keV; 
using compTB[2], we detected the presence of a blackbody 
component: the blackbody luminosity does not vary 
(0.3 x 10 38 erg sec -1 ) and its emission region is of a di- 
mension compatible with the NS and/or TL. 




10 20 
Energy (keV) 

Fig. 7. Unabsorbed EF(E) spectrum of GX 349+2 for HR<7. 
The dotted lines are the two separated compTB components: 
dash for compTB[l], dash-dot compTB[2]. 

Table 7. Best-fit parameters of the model wabs • (compTBfl]) 
for the source GX 3 + 1 . Errors are computed at 90% confidence 
level for a single parameter. 



parameter 


HR<8 


8<HR<9.5 


HR>9.5 




compTB[l] (thermal: log .A 


= 8,5 = 0) 


kT s a 
kT e a 
a 

T C 


n q+o-09 

n 7+0.1 
z - ' -0.1 

J-^-0.2 

4.4 


1 - u ^-0.06 
q+0.1 
°-0.2 
1 q+0-2 

3.2 


1 1Q+U.U7 
± - lo -0.06 

2.8l°i 
1.9 


CAF C 


1.4 


1.2 


1.1 


T b,c 
iJtot 


0.3 


0.2 


0.2 


X'fdof 


65/106 


116/114 


100/108 



a In keV. 

6 In units of 10 38 erg sec -1 , in the 0. 1-200 keV energy range. 
c Computed as reported in section |2~4l 



3.5. GX3+1 

In the case of GX 3+1, we identified three intervals of HR in 
which we recognized different phases of spectral evolution in 
the JEM-X1 and IBIS/ISGRI spectra: the three mean spectra 
obtained (HR<8, 8<HR<9.5 and HR>9.5) all represent the 
high/soft state. The spectral fitting was carried out on all spec- 
tra with a single component model, a thermal compTB (see 
Table|7] Fig.©: 

- as usual with compTB[l] we describe the TC process; sim- 
ilarly to previous sources, for increasing HR values, we see 
a decreasing of the TC efficiency (decrease of CAF and in- 
crease of a); 

- the seed photon temperature kT s is of the order 1 keV; 

- the plasma temperature kT e is of the order 3-4 keV. 



4. Discussion 

In this work we have studied the spectral evolution of five 
NS LMXBs using th e new Comptonization model compTB 
(Fari nelli et al]|2008l) . The compTB model is formed by two 
components: one component is a pure blackbody directly seen 
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Fig. 8. Unabsorbed EF(E) spectrum for the HR interval HR<8 
for GX 3 + 1 . The overall spectrum is well described by a single 
component, compTB[l], see text. 



by the observer while the other one takes into account the con- 
tribution of thermal and bulk Comptonization of the blackbody 
seed photons. To fit the INTEGRAL broad band spectra we ap- 
plied two compTB models. 

Every source has shown a different spectral evolution: for 
GX 5-1 and GX 13+1 increasing values of the HR correspond 
to spectral changes from the intermediate state (high-energy 
hard tail), to the high/soft state (no hard tail) while in the other 
three sources we did not observe dramatic spectral changes de- 
spite a progressive softening for increasing accretion rate (i.e. 
HR, see below). 

In all sources, regardless of HR, the total energetic bud- 
get is dominated by the thermal Comptonization component 
of kT s < 1 keV BB-like seed photons off warm electrons with 
temperature kT e typically of 3-4 ke V and optical depth of a few 
(t < 5). This component, described with a thermal compTB 
(equivalent to a compTT) is likely associated to the region be- 
tween the geometrically thin accretion disk and the NS surface, 
the so-called Transition Lay er, where most of the gra vitational 
energy release takes place ( Titarchuk & Fioritoll2004 , F08). 
When a second compTB model is required by the data, a sec- 
ond hotter seed photon population (~ 1-2 keV) can be iden- 
tified, interpreted as the BB-like contribution, from the inner- 
most region of the system, the NS and the inner TL (Fig. [TJ. 

A rough estimation of the apparent emitting radius of this 
BB region provides in fact results which are consistent with 
a NS and/or TL radius, with ^2-12 km. In the case of 
GX 5-1 for which we have the best statistics, we can see 
(Table [3} an increase of the emitting BB region with the HR 
(hence M, see next section). 

Where a single compTB model has been applied, we ob- 
serve that kT s has values very near to 1 keV, unlike the other 
cases where kT s is either below this energy (<0.6keV) or 
higher (always > IkeV). In the single compTB model case, 
it is likely that a part of the seed photons come from the disk 
(very soft photons <0.5 ke V) but emission from the TL and NS 
surface can contaminate the seed photon population increas- 
ing its overall temperature (~1 keV). As none of the observed 
sources reached the hard state (see, e.g. Fig. 2 in P06) where 
an increased height-scale of the TL may provide a high geo- 
metrical coverage of the central compact object, the lack of a 



second component in the case of a single compTB component 
in the fit is likely due to poor statistics in addition to small en- 
ergy coverage in the band where most of the BB-like emission 
is present. 

In terms of the observable quantities which may trace the 
source spectral evolution, it is important to focus on the a- 
parameter, related to the slope of the Green's function of the 
Comptonization region. As already pointed out in F08, for in- 
termediate state sources, we are dealing in fact with two ob- 
servable indexes, one related to the pure TC process, while the 
second one related to the TC+BC efficiency. The first index, 
only depends on the system temperature and optical depth and 
is actually directly related to the Comptonization parameter j0, 
while the second index contains an additional dependence on 
the bulk parameter 5 (see Eq.[23] in TMK97). In both cases, 
the slope of the Green's function depends on the first eigen- 
value of the spatial problem for a bounded medium. In fact, 
when BC is present, it is possible somewhat to treat the two re- 
gions (internal bulk-dominated TL, external TC dominated TL 
plus innermost accretion disk) in a separate way, each one with 
its own spatial configuration and associated eigenvalue prob- 
lem and derived spectral index. This translates into two differ- 
ent observable quantities (the thermal index a, plus the TC+BC 
index). For systems where BC is suppressed, or strongly re- 
duced, 5 goes to zero, and according to Eq.[13] in TMK97, the 
two indexes tend to be coincident and one single TC compo- 
nent is observed. 

The a parameter measures the slope of the spectrum 
and can quantify the effective efficiency of the overall 
Comptonization processes i.e. how much energy has been ex- 
changed between the photon field and the electrons. In the first 
compTB (the thermal one), a becomes higher as HR increases: 
we can always observe the spectral evolution in the direc- 
tion of a progressive softening, a smaller quantity of energy is 
gained by the photon field reducing the maximum energy of the 
spectra. Indeed, the energy gained by thermal Comptonization, 
CAF, decreases as a increases, giving a physical evaluation of 
the energy gain through the TC process. The a value related to 
the mixed compTB is generally higher with respect to the ther- 
mal one. It encloses both the BC and TC contributions coming 
from the TL, the NS itself and possibly the inner accretion disk. 

4. 1. The transient nature of the hard tail 

Bulk Comptonization acts in the TL, where the orbits of the in- 
falling matter deviate from the Keplerian behaviour to adjust to 
the NS rotation speed. Here the matter, during its bulk motion 
(vb ~ 0.1c, Titarchuk & Farinelli in prep.), interacts with the 
photon field that originates in the inner disk, in the TL itself 
and from the impact of matter onto the NS. 
Since the physical conditions of the interaction region strongly 
influence the BC, this produces the transient high energy hard 
tail behaviour. Indeed the local accretion rate of the source de- 



The 

fcT. 



Comptonization parameter is defined as 
max(r e ,r e 2 ) where the optical depth of the scattering 



V = 

region in NS LMXBs is typically t% » 1. 
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termines the resulting radiation pressure in the inner TL and in 
turn influences its hydrodynamical configuration. 

At very high values of accretion rate , the radiation pressure 
behaves like an intense photon wind (IBradshaw et al.1 12007) 
able to break efficiently the accreting matter, abating the t>&, 
hence suppressing (or strongly reducing) the Fermi first-order 
process in photon energy gain which is responsible for genera- 
tion of the hard X-ray tail. In the compTB model, the quantity 
that allows us to evaluate the hard tail- M link is the bulk pa- 
rameter 6 which is defined as 



(AE 



Bulk 



1 - I 



(5) 



where I = L/Lsdd and m = M / Msdd are local luminosity 
and accretion rate in Eddington's units, respectively, and = 
kT e /m e c 2 is the adimensional electron temperature! 
Increasing the accretion rate leads the system to a state closer 
to thermal equilibrium between electrons and the background 
radiation field, through a more efficient Compton cooling of the 
electrons themselves (lower average energy gain of the photons 
and CAF decreasing). 

In this work this spectral evolution has been observed for 
increasing values of HR with a slow spectral softening start- 
ing from the intermediate state (low HR) to the high/soft state 
(high HR values). Given the spectral evolution observed in all 
the sources, we can reasonably associate the increase of HR to 
an increasing accretion rate As a l ready pointed out by previous 
autho rs (e.g. iHoman et al.l 120071 : iLin et al. 2009; van der Klis 
200 lb . the role of M in the spectral evolution of these sources 
is not straightforward. Furthermore, from the observational re- 
sults of a detailed time-resolved BeppoSAX spectral analysis of 
the Z source Cyg X-2, F09 proposed to split the general defi- 
nition of the accretion rate M, into two quantities, Mdisk an d 
Mtl> obeying the condition A/tl Si A/disk- The reason for this 
resides in the fact that the total measured source bolometric lu- 
minosity does not follow the expected M trend from the source 
position in the hardness intensity diagram, and thus cannot be 
used as a tracer of the total mas s accre tion rate. Moreover, as 
also reported in IBradshaw et al.l ( 120071) . for high (Eddington- 
like) accretion rates, a strong radiation pressure may originate 
in the accretion disk, which may eventually eject part of the ac- 
creting material, so that the mass flow at the disk/TL radius is 
actually lower than the one flowing trough the Roche lobe. In 
this picture, the key role in the innermost system is played by 
Mtl- 

We can thus read the observed spectral variability of our 
sources as the interplay of TC and BC processes according to 
the aforementioned considerations and increasing HR: 

1. for a low Mtl, the hard tail results to be absent or too 
weak because of a lack of a sufficiently high number of BC 
scatterings able to elevate the hard tail contribution over the 
sensitivity instrument threshold. The observed spectrum is 
due mainly to the TC of soft seed photons from the disk and 



4 We note that the correct express i on o f the 5 parameter is the one 
written in the text. In lFarinelli e t al. (2008), Eq.[2] erroneously indi- 
cates a square root in the right term. 



NS in the hot corona plasma (tens of keV). This is the case 
of low-dim Atolls not considered here (see F08); 

2. for higher Mtl, the corona starts to be efficiently cooled 
and the overall spectrum is dominated by TC in a cold 
corona (kT e of a few keV). BC may not be strong enough 
to be detected (GX 3+1, GX 9+1, GX 349+2, but see also 
point 4); 

3. for increasing Mtl, we can observe the intermediate state 
(GX 5-1, GX 13+1 in this work), in which the hard tail 
is visible. Its importance and intensity are strongly influ- 
enced by the accretion rate value. Indeed only within a nar- 
row range of Mtl, the number of BC interactions is high 
enough to produce the hard tail and at the same time the ra- 
diation pressure is not too high to reduce the BC efficiency. 
The fine tuning of these two quantities could be the rea- 
son why we observe a little number of sources hosting NS 
displaying hard tails, and even more, only a given spectral 
state within a given source (the intermediate state, corre- 
sponding to the so-called Horizontal Branch of Z sources); 

4. in case of very high values of Mtl, the hard tail disappears 
or its contribution results to be negligible being under the 
sensitivity threshold. The pressure in the TL increases be- 
cause a large amount of matter releases a lot of energy at 
the impact onto the NS or within the TL itself and this can 
inhibit the bulk motion. The source moves to the high/soft 
spectral state in which the spectrum is mainly defined by 
the TC process plus eventually a variable contribution of a 
blackbody component emitted from the TL and/or the NS 
surface (GX 3+1, GX 9+1, GX 349+2). 

In the case of GX 5-1 for which we have the best statistics 
and HR evolution, we can relate the increasing trend of the 
Rbb value (size of the BB emitting region) to the increasing 
values of HR/M : at large values of local accretion rates, i.e. ra- 
diation pressure, bulk motion is inhibited and the matter falls 
at lower speed, following a complex behaviour for which it 
may arrive on the NS surfa ce spreading over a larger area 
dlnogamov & Sunvaev 1999b . The large amount matter cover- 
ing the NS can produce a layer able to suppress the radio emis- 
sion (P06). 

4.2. The absence of the hard tail in most bright Atoll 
"GX" sources 

From the spectral states observed and the information obtained 
by the fitting process, it is difficult to say whether the per- 
manent absence of hard X-ray tail in sources like GX 3+1, 
GX 9+1 is due to a locally high-enough accretion rate which 
is able to suppress BC or in fact to the opposite case, where BC 
is present but too low to be detected. In these sources such a 
transient feature has in fact never been observed, even with the 
past missions. 

In GX 5-1, a transient hard X-ray tail was instead pre- 
yiously detected b y Ginga dAsai et alJll994h and INTEGRAL 
JPaizis etalJ E005). in GX 13+1 it was detected for the first 
time by INTEGRAL (P0 6), while in GX 349+ 2, the feature was 
observed by BeppoSAX jDi Salvo et al.ll200ll) when the source 
was out-of-flare. Also taking into account previous results on 
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other Z sources (GX 17+2, Sco X-l, GX 340+0, see e.g. P06 
for a complete list of references) it is currently difficult to pro- 
vide a global self-consistent picture which may explain this 
variate phenomenology. We still do not have an observable pa- 
rameter which may help to unambiguously determine, at least 
phenomenologically, which is the threshold for triggering the 
hard X-ray tail. 

In the accretion scenario where the production of high- 
energy photons is inhibited by high local levels of the accretion 
rate (bulk stopping through radiation pressure), one should look 
with particular care at the contribution of the BB-like emission 
to the total luminosity. If this BB-like emission is due to en- 
ergy release of the accreting matter both to the NS surface and 
in the TL (through viscous dissipation), one would expect that 
hard tail quenching would be correlated to higher percentage 
of the BB contribution to the total luminosity. 

Using Lbb / L tot as a possible tracer of the local accretion 
rate, and thus of the BC efficiency, would have the advan- 
tage of being distance-independent, a not secondary issue when 
comparing the estimated luminosity among sources. Such a 
quantity may however be significantly different among sources 
which show a hard X-ray tail, ranging from ~2% (in GX 13+1, 
this work) to ^25% in Cyg X-2 (F08). Of course, a satisfac- 
tory energy coverage below 1 keV is of prime importance to 
give unbiased estimations of this contribution. In this sense the 
reported values of L to t and Lbb / -^tot m the present work must 
be treated very cautiously, as the energy band covered by JEM- 
XI starts from 5keV, thus above the threshold where most of 
BB energy is emitted. 

With this prescriptions in mind, following the spectral evo- 
lution of GX 5-1 with increasing HR, we note that the presence 
of the hard X-ray tail is anti-correlated with the total source 
luminosity and the percentage of the energetic budget carried- 
out by the BB-like component, even though the latter change is 
marginal (see Fig.@]and Tabled- The other sources for which 
we can detect two components in the spectrum along the ob- 
served HR range, but with no bulk contribution detected, un- 
like for GX 5-1, are GX 9+1 and GX 349+2 (see Figs.|6]and 
[7]). In this case case the BB-like component contributes about 
20% and 10% of the source luminosity, respectively. Little can 
be said though for GX 3+1 or GX 13+1 for which we do not 
have the double component evolution along with HR. 

5. Conclusions 

Despite an increasing amount of theoretical and observational 
results, the picture giving rise to the presence of variable PL- 
like hard X-ray emission in NS LMXBs is far from being 
completely clear. Investigating the physical conditions of the 
innermost region of these sources is very critical. The bulk 
Comptonization approach, together with the physical scenario 
proposed, can be considered an important step forward in try- 
ing to understand what is at the origin of the transient hard tail 
behaviour. Nevertheless, it is currently not easy to satisfactorily 
explain why sources like GX 9+1 and GX 5-1 that do have a 
similar X-ray spectrum and related parameters (kT e , kT s of the 
disk and NS-TL, r) yet differ for such an important feature as 
the presence of the transient hard tail. 



Even if we were able to quantify Mtl (disentangling it 
from the overall M) we would face the problem of understand- 
ing what makes Mtl change in the different cases i.e. what is 
the origin of the inflow anisotropy. In fact, a theoretical and/or 
observational quantity which may unambiguously trigger or 
dump the hard X-ray emission has not been identified yet. Of 
course it is possible that we are dealing with a multi-parametric 
problem. 

A possible highly speculative explanation could reside in 
the different composition of the disk: in Atoll sources the or- 
bital period is in general relatively short (^5 hr) and in the Z 
sources for which it is known it is long (^10 hr). This would 
indicate that Z sources may contain an evo lved companion 
star (e.g. Cyg X-2, GX 1+4 iLiu et all l2007l and references 
therein) w hile Atolls con tain a main sequence star (e.g. Aql X- 
1 , Ser X-l ILiu et al.ll2007l and references therein). This is con- 
sistent with the fact that the number of known Atoll sources 



(25) i s higher than the number of known Z sources (8. lLiu et al 



2007) since the main sequence phase of a star is much longer 



than the one o f an evolved star making the det ection of Atolls 
more likely dVerbunt & van den Heuvel 1995 ). The accretion 
disk coming from an evolved companion (Z-like system) would 
be Helium-rich (i.e. heavier) hence a higher luminosity would 
be required to stop the bulk flow, making the detection of BC 
more likely. 

The chemical composition of the disc could be inferred by 
the duration of type I X-ray bursts: long bursts are thought to 
be due to mixed H/He burning, triggered by thermally unsta- 
ble H e ignition, expected a t inferred near-Eddington accretion 
rates dKuulkers et alj|2007l) . However the observed behaviour 
is far from being interpreted within the current bursting the- 
ory due to the presence of unexpected short burst or even ab- 
sence thereof, in a regime where long bursts would be expected 



(Kuulkers et al. 2007). The reason of this could reside in the 



fact that in the case of high accretion rates, accretion is un- 
likely through a plain disc and the TL configuration may also 
have an effect on the bursting behaviour. 

Additional or concurring effects could reside in the pres- 
ence of a counter-rotating NS with respect to the Keplerian 
disk rotation, whose main effect could be the presence of a 
higher level of gravito-kinetic energy release in the TL, in the 
region where the angular velocity adjusts from the Keplerian 
regime to that of the slowly spinning NS. Such higher energy 
release would in turn produce a higher local radiation pressure 
gradient, more efficient in stopping bulk. Unfortunately, this is 
highly speculative and currently difficult to prove. A first ap- 
proach to tackle the problem would be a more detailed study of 
the spectral evolution together with the long-term behaviour of 
these sources (Savolainen et al, in prep.). 

Understanding the origin of the hard tails in NS LMXBs is 
important to make a decisive step forward, from phenomenol- 
ogy (power-laws) to physics. This step forward will not make 
us understand hard tails alone, on the contrary, hard tails could 
be the observational feature that gives us the means to un- 
ders tand the mechanism of ac creti on flows in g e neral, see 
e.g lShrader & Titarchukl d 19981) and lFerrigno et all d2009h for 
an application of bulk motion Comptonization in Black Hole 
LMXBs and High Mass X-Ray Binaries, respectively. 
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